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摘 要 考虑 详细 化 学 反应 及 输 运 过 程 ， 对 不 可 压 N-S 方程 直接 数值 求解 ,研究 了 低温 化 学 影响 下 的 火焰 水 力学 及 胞 状 不 稳定 
性 。 结 果 表 明 ， 预 热 后 已 发 生 低 温 化 学 反应 的 混合 气 ， 其 火焰 的 扰动 增长 率 明 最 大 于 未 发 生 低温 反应 的 混合 气 火 焰 。 通 过 对 预 
热 后 混合 气 火焰 的 关键 参数 分 析 ， 发 现 和 未 预 热 混 气 火焰 相 比 其 无 量 纲 释 热量 、Lewis 数 及 Zeldovich 数 都 显著 减 小 。 依 据 已 
理论 可 知 ， 在 Lewis 数 大 于 1 的 条 件 下 ， 上 述 参数 的 减 小 将 使 得 胞 状 不 稳定 性 对 扰动 增长 的 抑制 作用 减弱 。 因 此 ， 预 热 后 火焰 
的 增长 率 明 显 增 大 。 同 时 ， 在 低温 反应 中 ， 大 量 自由 基 和 其 它 小 分 子 的 形成 ， 使 得 预 热 后 混合 气 火焰 Lewis 数 减 小 ， 从 而 其 火 
焰 胞 状 结构 的 深度 和 高 度 都 变 小 。 
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Abstract Numerical simulation on the incompressible N-S equation was completed to study the low temperature 
chemistry on flame instabilities by using the detailed chemical kinetics and transport property. The results show that 
the flame disturbance growth rate of the mixed gas involving low temperature chemistry (LTC) is significantly 
greater than that without LTC. By the analysis of the key parameters of the preheated mixed gas flame, it is found 
that the dimensionless heat release, Lewis number and Zeldovich number are significantly reduced compared with 
that of the flame without preheating. On the basis of the existing theory that the reduction of the above parameters 
makes the cellular instability decreased inhibitory effect on the growth of disturbances as Lewis number greater than 
1. As a result, the growth rate of flame is obviously increased after preheating. At the same time, a large number of 
free radicals and other small molecules are formed in the low temperature reaction, which makes the flame Lewis 
number of the preheated mixture reduced, and the depth and height of the flame cellular structure become smaller. 
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0 前 言 

预 混 火 焰 水 力学 不 稳定 性 和 胞 状 不 稳定 性 通 
常 耦合 出 现 。 当 Lewis 数 小 于 1 时 , 后 者 强化 前 者 ， 
大 于 1 时 则 抑制 。 水 力学 不 稳定 性 在 各 个 波 数 下 都 
会 发 生 ， 而 在 考虑 火焰 厚度 以 及 拉 伸 效应 后 ， 其 在 
很 小 波 数 下 是 不 发 生 的 。 上 述 不 稳定 性 诱发 了 火焰 
面 的 裙 皱 ， 从 而 增 大 了 火焰 传播 速度 。 这 是 燃烧 室 
火焰 自满 流 化 ， 以 及 爆燃 转 爆 震 的 重要 机 制 。 在 当 
前 先进 发 动机 内 ， 广 泛 采 用 增 压 及 贫 预 混 燃烧 。 由 
于 压力 的 增 大 , 预 混 气体 的 点 火 延迟 时 间 显 著 减 小 。 
特别 是 对 于 航空 煤油 、 汽 油 及 柴油 等 大 分 子 燃料 混 
合 物 ， 其 存在 明显 的 低温 化 学 反应 。 因 而 ， 高 压 及 
高 预 热 温度 使 得 其 低温 反应 强烈 ， 同 时 低温 点 火 延 
退 时 间 不 断 减 小 。 当 火焰 低温 反应 时 间 尺 度 和 燃烧 
室 大 尺度 流动 时 间 尺 度 相似 时 ， 即 出 现 低温 化 学 影 
响 下 的 火焰 动力 学 过 程 四 。 此 时 ， 未 燃 混 气 并 非 初 
台 氧 化 剂 以 及 燃料 ， 而 是 低温 反应 路 径 中 的 中 间 产 
物 ， 以 及 剩余 的 燃料 和 氧化 剂 。 在 此 之 前 ， 关 于 水 
力学 及 胞 状 不 稳定 性 的 研究 都 是 基于 高 温 燃 烧 的 ， 
而 没有 学 者 考虑 过 低温 反应 下 火焰 面 传播 时 的 不 稳 
定性 机 理 。 而 在 实际 燃烧 室 中 , 当 低 温 反应 发 生 时 ， 
混合 气 的 化 学 及 热 性 质 都 发 生 了 显著 变化 ， 火 焰 伟 
播 不 稳定 机 制 也 一 定 会 发 生变 化 。 本 文 就 将 对 贫 燃 
预 混 二 甲 醚 /空气 , 低温 化 学 反应 影响 下 的 火焰 传播 
不 稳定 性 进行 研究 。 


1 物理 及 数学 模型 


本 文 考虑 的 是 二 维 ， 非 稳 态 反应 流 。 对 不 可 压 
NS 方程 采用 分 步 投 影 法 进行 直接 求解 。 方 程 离散 
保持 四 阶 精度 ， 非 稳 态 项 采用 C-N 格式 。 化 学 反应 
机 理 采 用 已 在 DNS 中 成 功 使 用 的 ， 去 除 反应 刚性 
的 简化 机 理 握 .此 机 理 包 含 30 个 组 分 和 175 步 反 应 。 
计算 过 程 中 ， 本 研究 耦合 Chemkin 软件 包 B1， 考 虑 
详细 的 化 学 反应 。 为 了 精确 计算 组 分 输 运 ， 同 时 耦 
合 Transport 软件 包 外 。 和 之 前 应 用 此 机 理 的 DNS 
研究 不 同 ， 本 文 考虑 详细 的 组 分 输 运 过 程 。 研 究 的 
物理 模型 为 二 维 平面 火焰 ， 左 侧 为 未 燃 混 气 进 口 ， 
右 侧 为 已 燃气 体 出 口 。 上 下 边界 采用 周期 边界 条 件 。 
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计算 初始 ， 将 采用 Chemkin 的 Premix 程序 计算 得 
到 的 层 流 火焰 结构 ， 设 置 在 二 维 计算 域 中 ， 并 设置 
火焰 位 置 为 Sin 函数 分 布 ， 扰 动 幅 值 为 一 个 火焰 厚 


度 。 
计算 区 域 长 度 为 60 倍 火焰 厚度 ， 高 度 为 一 个 
扰动 波长 。 网 格 分 辨 率 为 600x120， 时 间 步 长 为 


1x10*s。 为 了 验证 本 程序 的 准确 性 ， 将 此 程序 计算 
得 到 的 一 维 火 焰 结 构 和 Chemkin 中 Premix 程序 计 
算 值 进行 了 对 比 ， 如 图 1 所 示 。 从 图 中 可 以 看 到 ， 
当前 程序 以 及 网 格 分 辩 率 ， 能 够 精确 地 计算 层 流 火 


焰 。 


Heat Release Rate / erg/cm3-s 
Temperature /K 


T T 有 T T 
1.47 1.50 1.53 1.56 1.59 1.62 1.65 
x/cm 


图 1 DME/ 空 气 火焰 在 常 压 和 当量 比 为 0.5 时 的 温度 及 释 热 


Fig. 1 Profiles of temperature and heat release rate for DME/air 


flame under ambient pressure and equivalence ratio of 0.5 
2 结果 与 讨论 


对 于 大 分 子 燃料 , 在 负 温度 系数 (NTC) 区 域 ， 
存在 显著 的 低温 化 学 反应 过 程 。 本 文选 取 了 反应 机 
理 较 简单 的 二 甲 醚 、 同 样 具有 这 种 低温 反应 特性 ， 
并 且 是 常用 的 替代 燃料 。 其 化 学 简化 机 理 ， 已 经 和 
详细 机 理 在 不 同 压力 和 当量 比 下 的 火焰 传播 速度 、 
点 火 延迟 时 间 以 及 包括 点 火 和 灭火 拐点 的 “S” 曲 线 
进行 了 对 比 。 对 比 结果 表明 ， 此 具有 30 组 分 ，175 
步 反 应 的 简化 机 理 ， 可 以 很 好 的 模拟 低温 及 高 温 反 
应 影响 下 的 火焰 行为 。 图 2 是 本 文 给 出 的 DME/ 空 
气 在 当量 比 为 0.5， 压 力 分 别 为 0.1MPa 和 1.5MPa 
时 的 点 火 延迟 曲线 。 图 2 中 实 线 为 高 温 点 火 ， 而 点 
线 所 代表 的 是 低温 点 火 。 在 NTC 区 域 低温 反应 显 
著 ， 从 而 出 现 了 点 火 延 迟 时 间 更 短 的 低温 点 火 。 在 


当前 先进 发 动机 中 ， 随 着 压力 和 预 热 温度 的 提高 ， 
可 以 看 到 高 温 点 火 延迟 时 间 缩 短 的 同时 ， 其 高 压 下 
低温 点 火 延 迟 时 间 更 短 。 因 此 ， 在 混合 气 进入 高 温 
反应 区 燃烧 之 前 ， 就 已 经 发 生 了 低温 化 学 反应 ， 甚 
至 已 经 发 生 了 低温 点 火 。 本 文选 取 初 始 温度 为 
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到 在 0.1s 处 发 生 低温 点 火 。 可 以 看 到 低温 点 火 过 程 
中 ， 有 大 量 的 HO; 生成 。 这 是 因为 ， 在 低温 反应 过 
程 中 ，DME 首先 和 0 反应， 生成 CH3OCH (了 R) 

和 HO,, 即 RHHO;=R+HO,， 同 时 OH 也 会 夺取 RH 
中 的 互 并 也 生成 R， 即 RH+HOH=R+HO。 生 成 的 了 


600K， 当 量 比 为 0.5 的 DEM/ 空 气 混合 气 ， 对 其 二 
维 平面 火焰 进行 基于 详细 反应 和 输 运 的 直接 数值 模 
拟 。 基 于 计算 ， 研 究 其 低温 反应 对 火焰 水 力学 及 胞 
状 不 稳定 性 的 影响 。 图 3 即 为 其 释 热 速率 以 及 重要 
中 间 组 分 随时 间 的 变化 。 


1atm 高 温 点 火 


1000/T K’' 
图 2 当量 比 为 0.5, DME/ 空 气 在 常 压 及 15 个 大 气压 下 ， 
点 火 延 迟 时 间 随 初始 温度 的 变化 。 


Fig. 2 Ignition delay time as a function of initial 


temperature for DME/air flames at equivalence ratio 
of 0.5 under ambient pressure and elevated pressure of 
1.3 MPa. 
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均 质 (0D) DME/ 空 气 在 NTC 区 域 自 燃 过 程 释 
热 速率 以 及 重要 中 间 组 分 随时 间 的 变化 


Fig. 3 Profiles of heat release rate and medium 


图 


[50 


molecules versus time for homogeneous DME/air 


flames at NTC and autoignition 


由 于 所 取 工 况 属于 NTC 区 域 , 因此 可 以 明显 看 


与 02 反应 ， 形 成 重要 的 可 逆反 应 R+O2=RO:。 在 低 
温 条 件 下 ， 正 向 反应 占据 主导 ， 从 而 RO; 后 续 引 发 
的 链 分 支 反 应 将 引起 低温 点 火 。 所 以 ， 大 量 的 HO， 
产生 ， 说 明 同 时 有 大 量 的 R， 以 及 RO;, 的 产生 ， 从 
而 发 生 低温 点 火 。 在 低温 点 火 发 生 时 ， 同 时 有 大 量 
的 H20, 产生 。 这 是 由 于 HO, 大 量 产生 ， 反 应 
RH+HO2=R+H2O2， 以 及 HO2+HOs=H202+O,， 都 将 
生成 大 量 的 H20,。 此 组 分 的 不 断 积累 , 在 温度 升 高 
后 发 生 反 应 HIO=2OH， 从 而 引发 高 温 点 Si 因此 ， 
在 高 温 点 火 前 ， 可 以 看 到 大 量 的 H20; 迅速 减少 ， 
而 OH 在 高 温 点 火 处 迅速 上 升 。 

在 没有 考虑 低温 反应 时 ， 初 始 温度 为 600K 的 
DME/ 空 气 预 混 火 焰 传 播 过 程 如 图 4 所 示 。 


0.035 


8a 0.015 0.02 0.025 
X/cm 


图 4 初始 温度 为 600K， 当 量 比 0.5 的 DME/ 空 气 火焰 
面 随时 间 演 化 过 程 
Fig. 4 Evolution of flame surfaces for DME/air flames at 


initial temperature of 600K and equivalence ratio of 0.5. 


在 初始 时 刻 , 平面 火焰 位 置 登 加 Sin 函数 扰动 。 
区 域 高 度 为 一 个 波长 ， 此 波长 是 根据 火焰 Lewis 数 
确定 的 最 不 稳定 波 数 下 的 波长 。 火 焰 在 水 力学 不 稳 
定性 的 作用 下 ， 发 生 皱 折 。 在 50 个 单位 时 间 后 ， 火 
焰 基 本 稳定 不 再 变化 。 这 是 由 于 ， 基 于 混合 物 热 扩 
散 系数 和 燃料 的 质 扩 散 系 数 确 定 的 火焰 Lewis 数 为 
1.68。 胞 状 不 稳定 性 此 时 起 到 稳定 火焰 的 作用 ， 即 
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在 突出 位 置 火焰 速度 减 小 ， 而 在 凹陷 位 置 火 焰 速 度 


增 大 。 当 凹 火焰 变形 不 断 增 大 ， 即 拉 伸 率 增 大 的 同 a 
时 ， 当 地 火焰 速度 在 扩散 - 热 作用 下 会 增 大 ,最 终 形 有 
成 稳定 的 胞 状 结 构 。 而 凸 火焰 由 于 其 在 低 拉 伸 率 下 ， 
火焰 的 减速 和 火焰 前 流 场 匹配 ， 从 而 很 快 就 形成 稳 | 
定形 态 。 为 了 进一步 , 研究 火焰 初始 扰动 A0， 随 时 ee 


间 的 增长 率 ， 本 文 给 出 了 基于 Matalon 等 中 和 de 


Landaul | 等 的 理论 给 出 的 解 ,以 及 计算 得 到 的 结果 。 pb J i 
从 图 5 中 可 以 看 到 ， 本 计算 结果 和 Matalon 等 给 出 Em 
的 色散 关系 得 到 的 增长 率 基本 一 致 ,而 Laudau 等 的 图 6 当量 比 0.5 的 DME/ 空 气 混合 物 , 初始 温度 600K， 
理论 由 于 没有 考虑 胞 状 不 稳定 性 的 影响 ， 其 增长 率 预 热 0.4s 后 建立 的 火焰 面 随时 间 演 化 过 程 
明显 偏 大 。 这 也 再 次 指明 , 在 当前 Lewis 数 情况 下 ， Fig. 6 Evolution of flame surfaces for DME/air flames at 
火焰 胞 状 不 稳定 性 显著 的 抑制 由 水 力学 不 稳定 性 引 initial temperature of 600K and equivalence ratio of 0.5 after 
发 的 扰动 增 大 。 preheated time of 0.4s 
全] 一 Matalon and Matkowsky (SVF) 在 经 历 15 个 单位 时 间 后 ， 火 焰 扰 动 逐渐 不 再 
Computona esuts 放大 。 同 时 ， 胞 状 结构 的 深度 ， 及 高 度 都 明显 小 于 


未 预 热 混 气 火焰 。 对 预 热 混 气 火焰 和 未 预 热 混 气 火 


28 
3? 0 焰 扰 动 增长 进行 比较 发 现 ， 预 热 后 的 混 气 火焰 其 扰 
aa 三 15 动 增长 率 明显 变 大 ， 如 图 7。 
~ 1 为 了 说 明 其 原因 ， 本 文 对 预 热 前 后 重要 参数 i 
有 行 了 分 析 ， 如 表 1 所 示 。 
CO 对 于 Zeldovich 数 ， 预 热 后 明显 减 小 。 说 明 预 
ee t 热 启动 了 低温 反应 过 程 ， 使 得 预 热 后 混 气 火焰 整体 
一 图 5 扰动 增长 的 理论 与 计算 值 比较 活化 能 变 小 ， 反 应 区 相对 增 厚 。 同 时 ， 无 量 纲 放 热 
Fig. 5 Comparison of theory and calculation for 量 Q 的 减 小 ， 说明 由 于 低温 反应 ,特别 是 低温 点 火 
the disturbance growth 的 发 生 ， 将 部 分 化 学 能 已 转换 为 热能 释放 。 
yy | 2.0 
© 为 了 研究 低温 化 学 对 上 述 火 焰 不 稳定 性 的 影 1, 600k 400ms 
响 ， 本 文平 面 火焰 左 端 进口 为 初 温 600K， 预 热 0.4s 
后 的 混合 气 。 从 图 6 可 以 明显 看 出 ， 预 热 后 的 混合 0 
气 火焰 ， 其 演化 过 程 和 图 4 明显 不 同 。 
ts 
0.0 


图 7 预 热 混 气 火焰 及 未 预 热 混 气 火焰 扰动 增长 比较 
Fig. 7 Comparison ofthe disturbance growth ofpreheated 


flames and flames without preheating 


表 1 未 预 热 和 预 热 0.4s 后 ， 火 焰 参 数 对 比 
Table 1 Comparison of flame parameters for preheated 0.4s and 


flame without preheating 


Ze 0 Le Ze 
未 预 热 | 7.66 2.035 1.65 0.76 
预 热 0.4s | 4.67 1.025 1.29 0.67 


对 于 Lewis 数 ， 本 文采 用 文献 中 的 方法 ， 基 于 
马克 斯 坦 长 度 的 计算 ， 得 到 火焰 的 Lewis 数 。 表 中 
可 以 看 出 , 预 热 后 混 气 火焰 的 Lewis 数 显 著 降低 了 。 
这 是 由 于 , DME 分 子 预 热 后 , 低温 反应 消耗 了 大 量 
燃料 的 同时 ， 产 生 了 小 分 子 自 由 基 ， 以 及 H2、CO 
和 CH4 等 较 小 的 分 子 。Le。 是 临界 Lewis 数 ， 小 于 
此 值 时 ， 胞 状 不 稳定 性 将 促使 扰动 的 增 大 。 可 以 看 
到 ， 预 热 后 混合 气 火焰 需要 更 小 的 Lewis 数 ， 胞 状 
不 稳定 性 才能 有 利于 扰动 增长 。 

上 述 这 些 火焰 关键 参数 的 变化 ， 将 使 得 预 热 后 
火焰 本 质 不 稳定 性 发 生变 化 。 基 于 Landau 等 及 
Matalon 等 的 理论 分 析 ， 本 文 给 出 了 无 量 纲 放 热量 、 
Lewis 数 及 Zeldovich 数 与 火焰 水 力学 和 胞 状 不 稳定 
性 的 关系 ， 如 图 8、9 所 示 。co 和 ci 分 别 表示 水 力 
学 不 稳定 性 及 胞 状 不 稳定 性 对 总 的 扰动 增长 率 的 贡 
献 。 可 以 看 出 ， 当 放 热 量 增 大 时 ， 火 焰 的 水 力学 不 
稳定 性 增强 。 然 而 ， 由 于 Lewis 数 大 于 1， 胞 状 不 
稳定 性 随 着 放 热量 的 增 大 ， 其 抑制 扰动 增长 的 程度 
变 大 。 特 别 是 ， 随 着 Lewis 数 的 增 大 ， 这 种 抑制 趋 
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Fig. 8 Preheated flame growth rates versus the dimensionless 


heat release rates 


图 9 预 热 混 气 火焰 增长 率 随 Ze 
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Fig. 9 Preheated flame growth rates versus Zeldovich numbers 


3 结 论 


本 文 针对 平面 二 维 贫 预 混 二 甲 醚 /空气 火焰 , 进 


势 更 明显 。 同 时 ， 当 Zeldovich 数 ] 


定性 对 扰动 增长 的 抑制 程度 更 强 。 


赠 大 ， 胞 状 不 稳 


因此 ， 和 预 热 后 


混 气 火焰 相 比 ， 未 发 生 低温 反应 的 混合 气 由 于 具有 
更 大 的 Lewis 数 、 放 热量 Q 和 Zeldovich 数 ,在 Lewis 
数 大 于 1 的 条 件 下 胞 状 不 稳定 性 抑制 扰动 增长 的 能 
力 更 强 。 所 以 预 热 后 混 气 火焰 具有 更 大 的 增长 率 。 

同时 ， 由 于 Lewis 数 的 减 小 ， 其 胞 状 深度 及 高 度 都 
更 小 ， 这 也 和 之 前 基于 一 步 总 包 反 应 和 简单 输 运 关 
系 ， 针 对 不 同 Lewis 数 下 高 温 燃 烧火 焰 的 数值 模拟 


研究 结果 相同 。 


行 了 直接 求解 耦合 详细 化 学 反应 和 输 运 过 程 的 N-S 
方程 的 数值 模拟 工作 。 采 用 的 化 学 反应 机 理 ， 能 够 
准确 的 描述 二 甲 醚 低温 反应 特性 。 同 时 ， 考 虑 详细 
的 输 运 过 程 , 增加 了 模拟 扩散 - 热 引 起 的 胞 状 不 稳定 
性 的 准确 度 。 文 章 将 已 发 生 低温 点 火 的 混合 气 建立 
的 火焰 和 未 发 生 低温 反应 的 DME/ 空 气 混合 气 火 
焰 ， 其 扰动 演化 过 程 进行 了 对 比 。 结 果 表 明 ， 预 热 
后 混合 气 火焰 的 扰动 增长 率 显 著 增 大 。 为 了 对 这 一 
现象 进行 分 析 ， 本 文 提 取 了 基于 拉 伸 火焰 和 平面 火 
焰 综 合计 算得 到 的 Lewis 数 以 及 火焰 的 Zeldovich 
数 和 无 量 纲 放 热量 。 上 述 参 数 对 比 结果 表明 ， 在 预 
热 后 混合 气 建立 的 火焰 中 ， 这 些 参数 都 显著 减 小 。 
基于 已 有 的 理论 分 析 可 知 ， 这 些 量 的 减 小 表明 火焰 
胞 状 不 稳定 性 在 Lewis 数 大 于 1 的 条 件 下 ， 将 抑制 


初始 扰动 在 水 力学 不 稳定 机 制 下 的 增长 。 从 而 ， 揭 
示 了 低温 化 学 影响 下 ， 火 焰 水 力学 和 胞 状 不 稳定 性 
的 作用 规律 。 同 时 ， 由 于 Lewis 数 的 减 小 ， 可 以 看 
到 预 热 后 混合 气 建立 的 火焰 ， 其 胞 状 结构 深度 和 高 
度 都 明显 减 小 。 


参考 文献 


[1] 


[2] 


[3] 


[4] 


[5] 


[7] 


Won S H, Windom B, Jiang B, et al., The Role of Low 
Temperature Fuel Chemistry on Turbulent Flame [J]. 
Combustion and Flame, 2014, 161: 475-483 

Bhagatwala A, Luo Z, Lu T F, et al., Numerical and 
experimental investigation of turbulent DME jet flames []]. 
Proceedings of the Combustion Institute, 2015, 35(2): 
1157-1166 

Kee R J, Rupley F M and Miller J A, Chemkin II: a 
FORTRAN chemical kinetics package for the analysis of 
gas-phase chemical kinetics, SAND89-8009, Sandia 
National Laboratories 1989 

Kee R J, Dixon-lewis G Warnatz J, et al.,A Fortran 
Computer Code Package For The Evaluation Of Gas-Phase, 
Multicomponent Transport Properties, SAND86-8246, 
Sandia National Laboratories, 1986 

Matalon M, Matkowsky B J, Flames as gasdynamic 
discontinuities [J]. Journal of Fluid Mechanics, 1982, 124: 
239-59 

Landau L D, Lifshitz E M, Fluid Mechanics [MI]. Oxford: 
Pergamon, 1959 

Bouvet N, Halter F, Chauveau C, et al. On the effective 
Lewis number formulations for lean 
hydrogen/hydrocarbon/air mixtures [J]. International 


Journal of Hydrogen Energy, 2013, 38: 5946-5960 


